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ABSTRACT 

As part of several Hi synthesis mapping studies of merging galaxies, we have mapped the tidal 
gas in the three disk-disk merger systems Arp 157 (NGC 520), Arp 220, and Arp 299 (NGC 3690). 
These systems differ from the majority of the mergers mapped in H I, in that their stellar and gaseous 
tidal features do not coincide. In particular, they exhibit large stellar tidal features with little if any 
accompanying neutral gas and large gas-rich tidal features with little if any accompanying starlight. On 
smaller scales, there are striking anti-correlations where the gaseous and stellar tidal features appear to 
cross. We explore several possible causes for these differences, including dust obscuration, ram pressure 
stripping, and ionization effects. No single explanation can account for all of the observed differences. 
The fact that each of these systems shows evidence for a starburst driven superwind expanding in 
the direction of the most striking anti-correlations leads us to suggest that the superwind is primarily 
responsible for the observed differences, either by sweeping the features clear of gas via ram pressure, or 
by excavating a clear sightline towards the starburst and allowing UV photons to ionize regions of the 
tails. If this suggestion is correct, only systems hosting a galactic superwind and experiencing a high- 
inclination encounter geometry (such that tidal gas is lifted high above the starburst regions) should 
exhibit such extreme differences between their H I and optical tidal morphologies. 

Subject headings: galaxies: individual (Arp 220, NGC 3690, NGC 520) — galaxies: interactions — 
galaxies: ISM — galaxies: peculiar — galaxies: starburst — 



1. INTRODUCTION 

Nearly 30 years ago, Toomre & Toomre (1972) elegantly 
demonstrated that the tails and bridges emanating from 
many peculiar galaxies may arise kinematically from dy- 
namically cold disk material torn off of the outer regions 
of galaxies experiencing strong gravitational interactions. 
Early spectroscopic studies of gas within the tidal tails of 
merging galaxies provided observational support for this 
hypothesis by showing the tails to have the kinematics ex- 
pected for a gravitational origin (e.g. Stockton 1974a,b). 
H I mapping studies are particularly well suited to such 
studies, as the tidally ejected disk material is usually rich 
in neutral hydrogen and can be traced to very large dis- 
tances from the merging systems (e.g. van der Hulst 1979; 
Simkin et al. 1986; Appleton et al. 1981, 1987; Yun et al. 
1994). Once mapped, the tidal kinematics can be used 
either alone, to disentangle the approximate spin geome- 
try of the encounter (Stockton 1974a,b; Mihos et al. 1993; 
Hibbard & van Gorkom 1996, hereafter HvG96; Mihos & 
Bothun 1998), or in concert with detailed numerical mod- 
els, to constrain the full encounter geometry (e.g. Combes 
1978; Combes et al. 1988; Yun 1992, 1997; Hibbard & 



Mihos 1995; Gardiner & Noguchi 1996). 

However, not all systems can be easily explained by 
purely gravitational models such as those used by Toomre 
& Toomre. For example, gravitational forces by them- 
selves should not lead to differences between stellar and 
gaseous tidal components. Numerical models which in- 
clude hydrodynamical effects do predict a decoupling of 
the dissipative gaseous and non-dissipative stellar compo- 
nents (e.g. Noguchi 1988; Barnes & Hernquist 1991, 1996; 
Weil & Hernquist 1993; Mihos & Hernquist 1996; Apple- 
ton, Charmandaris & Struck 1996; Struck 1997), but only 
in the inner regions or along bridges where gas orbits may 
physically intersect (see e.g. Fig. 4 of Mihos & Hernquist 
1996). Decoupling of the gaseous and stellar components 
within the tidal tails is not expected. 

Nonetheless, differences between the optical and gaseous 
tidal morphologies have been observed. These differences 
can be subtle, with the peak optical and H I surface bright- 
nesses simply displaced by a few kpc within the tails (e.g. 
NGC 4747, Wevers et al. 1984; NGC 2782 Smith 1994; 
NGC 7714/4 Smith et al. 1997; Arp 295A, NGC 4676B, 
and NGC 520 Southern tail, Hibbard 1995, HvG96), or 
they can be extreme, with extensive H I tidal features ap- 
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1987, 1990 hereafter HAM90; Armus, Heckman, & Mi- 
ley 1990; Lehnert & Heckman 1996), and are thought to 
arise when the mechanical energy from massive stars and 
supernovae in the central starburst is sufficient to drive 
the dense interstellar medium outward along the minor 
axis (e.g. Chevalier & Clegg 1985; Joseph & Wright 1985; 
Suchkov et al. 1994). Often, such starbursts are power- 
ful enough to drive a freely expanding wind of hot plasma 
completely out of the galaxy ("blowout"; HAM90). 

In the following subsections we briefly discuss what is 
known about the dynamical state of each system, and 
describe the differences between the stellar and gaseous 
tidal morphologies. Throughout this paper distances and 
other physical properties are calculated assuming Hq = 
75 kms" 1 Mpc -1 . 

2.1. NGC 520 

NGC 520 (Arp 157, UGC 966) is an intermediate-stage 
merger, with the two progenitor nuclei separated by 40" 
(5.8 kpc for D= 30 Mpc, 1" = 145 pc) and embedded 
within a common luminous envelope (HvG96 and refer- 
ences therein; see Fig. |j]a). There is a bright optical tidal 
tail stretching 24 kpc to the southeast (henceforth referred 
to as the S Tail) which bends sharply eastward and con- 
nects onto a broad optical plume 3 . This plume contin- 
ues to the north and west for 60 kpc before it appears to 
connect onto extended light surrounding the dwarf galaxy 
UGC 957 (outside the region plotted in Fig. [l|ei; see Stock- 
ton & Bertola 1980). 

The primary nucleus (the easternmost nucleus in 
Fig. ^]a) possesses a massive (~ 5 x 10 9 M Q ) 1 kpc-scale 
rotating molecular gas disk (Sanders et al. 1988; Yun & 
Hibbard 1999). An Hi disk is kinematically centered on 
this molecular disk, and extends to a radius of ~ 20 kpc 
(labeled "Inner Disk" in Fig. |l|b). Beyond this there is an 
intermediate ring of H I with a mean radius of ~ 30 kpc 
(i.e., the material which contains the feature labeled "N 
Clump" in Fig. [I]b), and a nearly complete outer ring of 
H I with a mean radius of 60 kpc, which extends smoothly 
through the dwarf galaxy 4 UCG 957 (only partially seen in 
Fig.Pflb). There is a kinematic and morphological continu- 
ity between the molecular gas disk, the inner H I disk, and 
outer Hi ring (Yun & Hibbard 1999), which suggests that 
all of this material is associated with the primary nucleus. 

The observations suggests that the NGC 520 interac- 
tion involved a prograde-retrograde or prograde-polar spin 
geometry: The linear morphology of the optical tail-to- 
plume system is typical of features produced by a disk 
experiencing a prograde encounter (i.e., the disk rotates 
in the same direction as the merging systems orbit each 
other). The disk- like morphology and rotational kinemat- 
ics of the large-scale HI and the lack of any aligned lin- 
ear tidal features, on the other hand, are more typical of 
polar or retrograde encounter geometries (i.e., disk rota- 
tion either perpendicular to or opposite the direction of 
orbital motion). Such encounters fail to raise significant 
tails (Toomre & Toomre 1972; Barnes 1988), and much of 

3 We follow the naming convention of Schombert et al. 1990 and refer to tidal features with flat intensity profiles as plumes, and ones with 
Gaussian profiles as tails 

4 The HI kinematics show that UGC 957 is kinematically associated with this outer gas ring, although it may lie slightly above or below it. 
Rudimentary numerical modeling has shown that it is unlikely that the UCG galaxy is responsible for the main optical features of the NGC 520 
system (Stanford & Balcells 1991). It is unclear whether this system is an interloper or was recently assembled from the surrounding ring of 
gas. 



parently decoupled from, or even anti-correlated with, the 
optical tidal features. It is this latter category of objects 
that we wish to address in this paper. In particular, we 
address the morphology of the tidal gas and starlight in 
the merging systems NGC 520 (Arp 157), Arp 220, and 
Arp 299 (NGC 3690). 

The three systems were observed as part of our on-going 
studies on the tidal morphologies of optically and IR se- 
lected mergers (Hibbard 1995, HvG96, Hibbard & Yun 
1996 and in prep.). These studies involve moderate res- 
olution [Ofwhm ~ 15") VLA Hi spectral-line mapping 
observations and deep optical B and R broad-band imag- 
ing with large format CCDs using the KPNO 0.9m (NGC 
520) and the University of Hawaii 88" telescopes. The 
Hi and optical observations, reduction, and data prod- 
ucts have been presented in Hibbard (1995) and HvG96 
for NGC 520, in Hibbard & Yun (1999, hereafter HY99) 
for Arp 299, and in Yun & Hibbard (2000; see also Hib- 
bard & Yun 1996) for Arp 220. We refer the reader to 
these papers for details of the observations and data re- 
duction. These systems are extremely disturbed, and we 
cannot hope to offer a full description of their peculiari- 
ties here. For more information we refer the reader to the 
above references. 



2. OBSERVED STELLAR AND GASEOUS TIDAL 
MORPHOLOGIES 

Figures |l]-^| show the optical and atomic gas morpholo- 
gies of each of the three systems discussed here. For 
NGC 520 and Arp 220 only the inner regions are shown 
in order to highlight the differences we wish to address. 
Panel (a) presents a greyscale representation of the op- 
tical morphology of each system with features of interest 
labeled. Panel (b) shows the H I distribution. Contours 
indicate the distribution of H I mapped at low-resolution 
{Ofwhm ~ 30"), whereas the greyscales show the Hi 
mapped at higher resolution (Ofwhm ~ 15"). The former 
is sensitive to diffuse low column density (Nui) neutral 
hydrogen, while the latter delineates the distribution of 
the higher column density H I. The central region of each 
H I map appears to have a hole (indicated by the dotted 
contours), which is due to H I absorption against the radio 
continuum associated with the disk- wide starbursts taking 
place in each galaxy (see Condon et al. 1990). In panel (c), 
we again present the optical morphology in greyscales, and 
the higher resolution H I distribution as contours. Finally, 
panel (d) presents a smoothed, star-subtracted i?-band 
image contoured upon a greyscale representation of the 
high-resolution H I map. 

In the final panels of Figs. [l]-|^ dashed lines labeled 
"Slice" indicate the locations from which H I and optical 
intensity profiles have been extracted; these profiles are 
plotted in Figure ||. Arrows labeled "Superwind" indicate 
the position angle (P.A.) of Ha or soft x-ray plumes, be- 
lieved to arise from a starburst-driven outflow or galactic 
superwind in each system. Such outflows are common in 
other IR bright starbursts (e.g. Heckman, Armus & Miley 
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Fig. 1. — HI and Optical Morphology of NGC 520. Panel (a) (top left) displays an i?-band image with features of interest labeled. 
Panel (b) (upper right) presents a greyscale and contour representation of the H I data. The contours indicate the distribution of 
H I mapped at low-resolution (Ofwhm ~ 25"), whereas the greyscales show the H I mapped at higher resolution (Ofwhm ~ 17"). 
The lowest contour is drawn at a column density of 5x 10 19 cm -2 , with successive contours a factor of 2 higher. The greyscales range 
from lxl0 20 cm -2 (white) to 4xl0 20 cm -2 (black). Panel (c) (lower left) presents the optical image in greyscales with contours 
from the high- resolution Hi data superimposed. The lowest contour is drawn at a column density of lxl0 20 cm -2 , with successive 
contours a factor of 2 higher. Panel (d) (lower right) shows contours from a star-subtracted i?-band image (lowest contour of 
27 mag arcsec -2 and contour interval of 1 mag arcsec -2 ) upon greyscales of the high-resolution Hi data ffrom lxl0 20 cm -2 to 



4x10 cm ). The dashed line labeled "Slice" indicates the locations of the intensity profiles plotted in Fig. [4p,. The arrow labeled 
"Superwind" indicates the direction (but not the extent) of the putative expanding superwind. 



the disk material remains close to its original rotational 
plane. 

Neither the intermediate nor the outer H I ring has 
an optical counterpart (fiR> 27 mag arcsec -2 ). Despite 
smooth rotational kinematics, this outer H I has a very 
clumpy and irregular morphology, with notable gaps near 



the optical minor axis (labeled "NE gap" and "SW gap" 
in Fig. |l|c). This figure shows that the outer Hi and op- 
tical structures are anti-correlated, with the peak H I col- 
umn densities (associated with the N clump) located to 
one side of the optical plume. In Fig. |l|d the H I clump 
appears to be bounded on three sides by the optical con- 
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tours. In Fig. ^a we present an intensity profile at the 
location indicated by the dotted line in Fig. [l]d, showing 
that the gas column density increases precisely where the 
optical light decreases. 

While the H I features exhibit a clear rotational kine- 
matic signature, the well-defined edges and nearly-linear 
structure of the optical plume suggests that its constituent 
stars are moving predominantly along the plume, rather 
than in the plane of the sky: any substantial differential ro- 
tation would increase the width of the plume and result in 
a more disk-like morphology. We therefore conclude that 
the gas rings and optical plume are both morphologically 
and kinematically distinct entities. This suggests that the 
observed gas/star anti-correlation is either transient (and 
fortuitous) or actively maintained by some process. 

A deep Ha image of NGC 520 shows plumes of ionized 
gas emerging both north and south along the minor axis 
and reaching a projected height of 3 kpc from the nucleus 
(HvG96). It has been suggested that this plume represents 
a star burst-driven outflow of ionized gas (HvG96, Norman 
et al. 1996). The position angle of this plume is indicated 
by an arrow in Fig. ||d (P. A. = 25?). This direction corre- 
sponds to the most dramatic Hi/optical anti-correlations 
mentioned above, and in the following we suggest that this 
region of the optical tail actually lies directly in the path 
of the out-flowing wind. 

2.2. Arp 299 

Arp 299 (NGC 3690/IC 694, UGC 6471/2, Mrk 171, 
VV 118) is also an intermediate stage merger, with two 
disk systems (IC 694 to the east, NGC 3690 to the west 
- see Figure Ha) in close contact but with their respec- 
tive nuclei separated by 20" (4.7 kpc for D= 48 Mpc, 1" 
= 233 pc). A long, narrow, faint (/Jr^ 26mag arcsec~ 2 ) 
tidal tail stretches to the north to a radius of ~ 125 kpc. 
H I imaging of this system by HY99 (see also Nordgren et 
al. 1997) shows a rotating gas-rich disk within the inner re- 
gions, and a pair of parallel H I filaments extending to the 
north. From the Hi morphology and kinematics, HY99 
deduce that Arp 299 is the result of a prograde-retrograde 
or prograde-polar encounter between two late-type spirals, 
with the inner H I disk associated with the retrograde disk 
of IC 694, and the northern optical tail and tidal H I fila- 
ments ejected by the prograde disk of NGC 3690. 

The parallel-filament or bifurcated morphology of the 
tidal H I is quite unlike that of the optical tail. The inner 
H I filament (so labeled in Fig. |2^>) is of lower character- 
istic column density (Njji^ 8xl0 19 cm~ 2 ) and is associ- 
ated with the low surface brightness stellar tail (Fig. |c). 
The gas in this filament has a more irregular morphology 
than that in the outer filament (e.g., the "gap" and "knot" 
in Fig. ||b), and much of this material is detectable only 
after a substantial smoothing of the data (Fig. |2|b). The 
outer filament is characterized by a higher H I column den- 
sity (iVffi~1.5x 10 20 cm -2 ) but has no optical counterpart 
(hr> 27.5 mag arcsec" 2 ). This filament is displaced by 
approximately 20 kpc (in projection) to the west of the 
inner filament for much its length, after which the fila- 
ments merge together in a feature labeled the "N Clump" 
in Fig. ^3. The parallel filaments have nearly identical 

5 To produce filaments as well separated as found in Arp 299 HY99 
warp in the progenitor disk. 



kinematics along their entire lengths, and join smoothly 
at the N Clump. This implies that these features form a 
single physical structure. 

Based on preliminary numerical simulations, HY99 sug- 
gest that a bifurcated morphology can arise quite natu- 
rally during tail formation. This occurs when the optically 
faint, gas-rich outer regions of the progenitor disk are pro- 
jected adjacent to optically brighter regions coming from 
smaller initial radii 5 (see also Mihos, 2000). However, this 
scenario does not explain why the inner filament, presum- 
ably drawn from optically bright but still gas-rich material 
within the optical disk of the progenitor, should lack ac- 
companying H I. 

As in NGC 520, there is an anti-correlation between the 
optical and gaseous column densities across the N clump, 
with the highest gas column densities (2-3x 10 20 cm~ 2 ) lo- 
cated on either side of the optical tail. This is illustrated 
in Fig. |4jo, where we plot a profile along the position in- 
dicated by the dotted line in Fig. ||d. The optical tail 
emerges above the N clump, and appears to curve exactly 
around the northern edge of the N clump (labeled "hook" 
in Fig. ||c;). Also labeled in Fig. are the three regions 
with anomalously high H I velocity dispersions noted by 
HY99 {a hi ~ 13-20 km s _1 compared with a m ~ 7-10 
km s _1 for the remainder of the tail; see Fig. 7d of HY99) ; 
we will will refer to these regions in the discussion (§ 3.4.2). 

Within the main body of Arp 299, vigorous star forma- 
tion is taking place, with an inferred star formation rate 
(SFR) of 50 Mgyr" 1 (HAM90). Recent X-ray observa- 
tions reported by Heckman et al. (1999) show evidence for 
hot gas emerging from the inner regions and reaching 25 
kpc to the north, which the authors interpret as evidence 
for a hot, expanding superwind. The position angle of this 
feature (P.A.=25?) is indicated by the arrow in Fig. ||d, 
and points towards the inner tidal filament and N clump. 

2.3. Arp 220 

Arp 220 (UGC 9913, IC 4453/4) is the prototypical ul- 
traluminous infrared galaxy with Lg— ic%im = 1-5 x W 12 Lq 
(Soifer et al. 1984). It is an advanced merger system with 
two radio and infrared nuclei separated by 0.9" (345 pc 
for D= 79 Mpc), and a bright optical plume extending 35 
kpc to the NW (Fig. |3^i). Each of the two nuclei has its 
own compact molecular disk. The two nuclear disks are 
in turn embedded in one larger 1 kpc scale molecular gas 
disk (see Scoville, Yun, & Bryant 1997, Downes & Solomon 
1998, Sakamoto et al. 1999 and references therein). The 
spin axis of the eastern nucleus is aligned with that of the 
kpc-scale disk while the western nucleus rotates in the op- 
posite direction. These observations suggest that Arp 220 
is the product of a prograde-retrograde merger of two gas 
rich spiral galaxies (Scoville, Yun, & Bryant 1997). 

An irregular disk-like distribution of neutral hydrogen 
extends over a 100 kpc diameter region surrounding the 
optical galaxy (Yun & Hibbard 1999a). The overall Hi 
kinematics indicates that this material has a component 
of rotation in the same sense as that for the eastern nucleus 
and the molecular gas disk, and opposite the rotation of 
the western nucleus. This suggests that the H I disk and 
eastern nucleus originated from the retrograde progeni- 
tor, while the the western nucleus and NW optical plume 

suggest that the bifurcation is exacerbated by a pre-exisiting gaseous 
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Fig. 2. — HI and Optical Morphology of Arp 299. Panel (a) (top left) displays an _R-band image with features of interest 
labeled. Panel (b) (upper right) presents a greyscale and contour representation of the H I data. The contours indicate the 
distribution of H I mapped at low-resolution (Ofwhm ~ 35"), whereas the greyscales show the Hi mapped at higher resolution 
(Gfwhm ~ 20"). The lowest contour is drawn at a column density of 2.5x 10 19 cm -2 , with successive contours a factor of 2 
higher. The greyscales range from lxl0 20 cm~ 2 (white) to 4xl0 20 cm -2 (black). Panel (c) (lower left) presents the optical image 
in greyscales with contours from the high-resolution H I data superimposed. The lowest contour is drawn at a column density of 
5xl0 19 cm -2 , with successive contours a factor of 2 higher. Panel (d) (lower right) shows contours from a star-subtracted _R-band 
image (lowest contour of 27 mag arcsec -2 and contour interval of 1 mag arcsec -2 ) upon greyscales of the high-resolution H I data 
(from lxl0 2O cm -2 to 4xl0 20 cm -2 ). The dashed line labeled "Slice" indicates the locations of the intensity profiles plotted in 
Fig. Wb. The arrow labeled "Superwind" indicates the direction (but not the extent) of the putative expanding superwind. 



(Fig. [|a) arose from the prograde progenitor. 

Because of the vigorous star formation occurring within 
Arp 220 (SFR=340 M yr -1 , HAM90), much of the Hi 
within the optical body of the system is seen only in ab- 



sorption against the bright radio continuum emission from 
the central starburst. Beyond this, the Hi has high col- 
umn densities (JV.h-j~1.5x 10 20 cm~ 2 ), but only to the NE 
and SW. Most notably, there are local H I minima to the 
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Fig. 3. — HI and Optical Morphology of Arp 220. Panel (a) (top left) displays an i?-band image with features of interest labeled. 
Panel (b) (upper right) presents a greyscale and contour representation of the H I data. The contours indicate the distribution of 
H I mapped at low-resolution (Ofwhm ~ 30"), whereas the greyscales show the H I mapped at higher resolution (Ofwhm ~ 18"). 
The lowest contour is drawn at a column density of 5xl0 19 cm -2 , with successive contours a factor of 2 higher. The greyscales range 
from lxl0 20 cm~ 2 (white) to 4xl0 20 cm~ 2 (black). Panel (c) (lower left) presents the optical image in greyscales with contours 
from the high-resolution Hi data superimposed. The lowest contour is drawn at a column density of 5xl0 19 cm -2 , with successive 
contours a factor of 2 higher. Panel (d) (lower right) shows contours from a star-subtracted i?-band image (lowest contour of 
26 mag arcsec -2 and contour interval of 1 mag arcsec -2 ) upon greyscales of the high-resolution Hi data (from lxl0 20 cm -2 to 
4xl0 20 cm -2 ). The dashed line labeled "Slice" indicates the locations of the intensity profiles plotted in Fig. tUc. The arrow labeled 
"Superwind" indicates the direction (but not the extent) of the putative expanding superwind. 



NW and SE (see gaps in Fig. ||b) . Comparison of the H I 
map with the optical image (Fig. ||c) shows that the NW 
gap occurs exactly at the location of the optical tail. The 
relationship between the optical and H I surface brightness 
levels across this feature are illustrated by an intensity pro- 
file measured along the dotted line shown in Fig. [||d, and 
plotted in Fig. |c. As in NGC 520 and Arp 299, the gas 
column density increases precisely where the optical light 



from the tail begins to fall off. There is a similar H I gap 
to the SE, but in this case there is no corresponding op- 
tical feature associated with it. At even larger radii, the 
H I is more diffuse (iVjj/~3x 10 19 cm~ 2 ) and has no optical 
counterpart down to [Ir=27 mag arcsec -2 . 

An X-ray image obtained with the ROSAT HRI camera 
(Heckman et al. 1996) reveals an extended central source 
that is elongated along P.A.=135? (indicated by arrows in 
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Fig. ||). A deep Ha + [Nil] image of Arp 220 reveals ion- 
ized gas with a bright linear morphology at this same po- 
sition angle (Heckman, Armus & Miley 1987). The optical 
emission line kinematics are suggestive of a bipolar outflow 
(HAM90), and the physical properties of the warm and 
hot gas strongly support the superwind scenario for this 
emission (HAM90, Heckman et al. 1996). As in NGC 520 
and Arp 299, the position angle of the putative expand- 
ing superwind is in the same direction as the H I minima, 
i.e. NW and SE. 

3. DISCUSSION 

Figures |l]-|] provide evidence for both small- and large- 
scale differences in the distributions of the tidal gas and 
stars in these three systems. The small-scale differences 
are of the type illustrated in Fig. [|, whereby the gas col- 
umn density falls off just as the optical surface brightness 
increases at various edges of the tidal features. In NGC 
520 and Arp 220, the large-scale differences are between 
the outer H I rings and disks (which have no associated 
starlight) and the optical tails and plumes (which have no 
associated H i) . Although these features are kinematically 
decoupled at present (with the gas rings and disks pre- 
dominantly in rotation and the optical tails and plumes 
predominantly in expansion), it is possible that they had 
a common origin and have subsequently decoupled and 
evolved separately. In Arp 299, on the other hand, the H I 
filaments and optical tail have similar morphologies and 
continuous kinematics and are therefore part of the same 
kinematic structure. In this system we believe the bifur- 
cated tidal morpholog y re sults from a progenitor with a 
warped gaseous disk (§2.2 & HY99), and we seek to un- 
derstand why the inner filament is gas-poor, given that its 
progenitor was obviously gas-rich. 

In this section we investigate a number of possible ex- 
planations for these observations. In particular, we discuss 
the possible role played by: differen ces in the initial radial 
distribution of the gas and stars (§|3.lD, dust obscuration 
(§ |3.2| ), kinematic decoupling of the gas due to collisions 
within the developing tidal tail ( §|3.3| ), ram pressure strip- 
ping of th e gas, either by a halo or by a galactic scale 
wind (§3.4), and photoionization of the gas, either by the 
starburst or by local sources (§^). 

3.1. Differences in the Radial Distribution of Gas and 

Starlight 

In interacting systems the H I is often more widely dis- 
tributed than the optical light (see, e.g. the H I map of the 
M81 system by Yun et al. 1994; see also van der Hulst 1979; 
Appleton, Davies & Stephenson 1981). These gas-rich ex- 
tensions frequently have no associated starlight down to 
very faint limits (e.g. Simkin et al. 1986; HvG96). A natu- 
ral explanation is that such features arise from the H i-rich 
but optically faint outer radii of the progenitor disks. The 
relatively short lifetimes of luminous stars and the larger 
velocity dispersions of less luminous stars, especially with 
respect to the gas, will further dilute the luminous content 
of this material, and the H i-to-light ratio of the result- 
ing tidal features will increase with time (Hibbard et al. 
1994). Gaseous tidal extensions with very little detectable 
starlight would seem to be the natural consequence. The 
outer H I rings in NGC 520 and Arp 220 and the gas-rich 
outer filament in Arp 299 are all likely to have arisen in 



this manner. 

However, gas-rich outer disks cannot give rise to gas- 
poor optical structures, such as the optical plume in NGC 
520, the optical tail in Arp 220, or the inner filament in 
Arp 299. Since these features presumably arise from opti- 
cally brighter regions of the progenitor disks (regions which 
are characterized by H I column densities higher than that 
of the outer disks) one would have expected a priori that 
these features should also be gas-rich. It is possible that 
the disks which gave rise to the plumes in NGC 520 and 
Arp 220 were gas-poor at all radii. However, this would 
not account for the discontinuities in the outer gaseous 
features that project near these optical features (i.e., NE 
gap in NGC 520 and NW gap in Arp 220). We therefore 
seek other explanations for these structures. 

3.2. Effects of Dust Obscuration 

The correspondence between rising gas column density 
and falling optical surface brightness (Fig. suggests that 
dust associated with the cold gas may attenuate the optical 
light. To address this possibility, we calculate the expected 
extinction in the R— band for a given column density of 
Hi. We adopt the Milky Way dust-to-gas ratio deter- 
mined by Bohlin, Savage, & Drake (1978; N HI /E(B-V)= 
4.8x 10 21 cm -2 mag -1 ), which is supported by direct imag- 
ing of the cold dust in the outer regions of eight disk galax- 
ies (Alton et al. 1998). This is combined with the Galactic 
extinction law of O'Donnell (1994; A R /E(B - V) = 2.673, 
from Table 6 of Schlegel et al. 1998) to yield an expected 
extinction in the i?-band of Ar — 1 8xU ^f/ cm ~s mag. 

From Fig. [|, the peak H I column densities on either 
side of the optical features are ~3xl0 20 cm -2 . The pre- 
dicted extinction is therefore of order 0.2 mag in the R- 
band. From Fig. |i| we see that the mean light level drops 
by about 1.0 mag arcsec -2 for Arp 299 (from 26.5 mag 
arcsec -2 to below 27.5 mag arcsec -2 ), about 1.5 mag 
arcsec -2 for NGC 520 (from 25 mag arcsec -2 to below 26.5 
mag arcsec -2 ), and by about 2.5 mag arcsec -2 for Arp 220 
(from 23.5 mag arcsec -2 to below 26 mag arcsec -2 ) along 
the extracted slices. To produce this amount of extinction, 
the tidal gas would have to have a dust-to-gas ratio that 
is ten times that in the Milky Way. 

The above analysis assumes that the measured neutral 
gas column density represents the total gas column density. 
However, the sharp drop in H I column density observed 
in many tidal features (HvG96, Hibbard & Yun in prepa- 
ration) suggests that the tidal gas may be highly io nize d 
by the intergalactic UV field (see also references in §3.5). 
Since large dust grains should survive in the presence of 
this ionizing radiation, the opacity per atom of neutral hy- 
drogen (Ar/Nhi) should increase in regions of increasing 
ionization fraction. Observations of NGC 5018 (Hilker & 
Kisslcr-Patig 1996), in which blue globular clusters are ab- 
sent in a region underlying an associated H I tidal stream, 
may support a high A^/Nhi ratio for tidal gas. Never- 
theless, the lack of obvious reddening of the B — R colors 
along the slices in Arp 299 and NGC 520 (Hibbard 1995; 
HY99) argues against a much higher extinction in these 
regions. 

We conclude that extinction might be important for 
shaping the morphology of the faintest optical features 
(e.g., the "Hook" and the end of the optical tail of Arp 




Fig. 4. — Intensity profiles of i?-band surface brightness (uji) and H I column density (beam-averaged values measured from the 
low resolution data) taken along the dashed lines in Figs. |j- |3| (a) (top) NGC 520, across the northern plume, (b) (middle) 
Arp 299, across the tip of the northern tail, (c) (bottom) Arp 220, across the NW plume. Gas column densities are measured in 
Mq pc~ 2 (where 1 Mq pc -2 = 1.25xl0 20 cm~ ), and optical surface brightnesses in mag arcsec -2 . 



299, Fig. H which has /j,r near the detection limit of 28 
mag arcsec -2 ), but is insufficient to greatly affect the over- 
all tidal morphology. However, an anomalously high tidal 
dust-to-gas ratio remains a possibility. This question could 
be resolved by the direct detection of cold dust in tidal tails 
with sub-millimeter imaging. 

3.3. Collisions within Developing Tidal Tails 

During the tail formation process, the leading-edge of 
the tail is decelerated with respect to the center of mass 
of the progenitors, while the trailing-edge is accelerated, 
and the two edges move towards each other (see Toomre 
& Toomre 1972, Fig. 3). Eventually, the two edges appear 
to cross, forming a caustic (Wallin 1990; Struck-Marcell 
1990). In most cases, the caustics are simply due to pro- 



jection effects. Only for low-inclination encounters will 
these crossings correspond to physical density enhance- 
ments, and numerical experiments suggest that in these 
cases the density will increase by factors of a few (Wallin 
1990). It has been suggested that collisions experienced 
by the crossing tidal streams in such low-inclination en- 
counters may lead to a separation between the dissipa- 
tional (gas) and non-dissipational (stellar) tidal compo- 
nents (Wevers 1984; Smith et al. 1997). 

The present data do not allow us to directly address 
this question, since the kinematic decoupling presumably 
took place long ago. However, several arguments lead us 
to suspect that this collisional process is not important 
in tidal tails: (1) large scale decoupling between the stel- 
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lar and gaseous tidal morphologies is not seen in many 
systems known to have experienced low-inclination en- 
counters (e.g. NGC 4038/9, "The Antennae", Hibbard 
et al. in preparation; NGC 7252 "Atoms for Peace" Hib- 
bard et al. 1994; NGC 4676 "The Mice" HvG96); (2) the 
broad plume-like morphologies of the optical features in 
Arp 22 an d NGC 520 suggest rather inclined encounters 
(§§2.1, 2.3), in which case wide-spread collisions are not 
expected; and (3) the parallel filaments in the Arp 299 
tail have identical kinematics, whereas one would expect 
kinematic differences between the stripped and unstripped 
material. 

Therefore while gaseous collisions and dissipation might 
result in differences between gas and stars during tidal de- 
velopment (particularly along tidal bridges, where the gas 
streamlines are converging; e.g. Struck 1997; NGC 7714/5 
Smith et al. 1997; Arp 295 HvG96), we believe that they 
are not likely to lead to a wide-spread decoupling in the 
outer regions. 

3.4. Ram Pressure Stripping 

If the tidal features pass through a diffuse warm or hot 
medium, or if such a medium passes through the tidal fea- 
tures, it is possible that the tidal gas exchanges energy 
and momentum with this medium due to collisions. Such 
effects have been proposed to explain the stripping of the 
cool interstellar medium from spiral galaxies as they move 
through the hot IGM in clusters (Gunn & Gott 1972), 
and is referred to as Ram Pressure Stripping (RPS). Tidal 
features should be relatively easily stripped, as they lack 
the natural restoring forces present in disk galaxies, except 
possibly at a small number of self-gravitating regions. In 
this case, the momentum imparted due to ram pressure 
is simply added to or subtracted from the momentum of 
the gaseous tidal features, and a separation of stellar and 
gaseous components might be expected. 

In the next two subsections, we investigate two possi- 
ble sources for ram pre ssure : an extended halo associated 
with the progenitor s (§[3. 4.1 ); and an expanding starburst 
driven superwind (§3.4.2). 



3.4.1. RPS from Extended Halo Gas 

Our own galaxy is known to have an extended halo of 
hot gas (Pietz et al. 1998). The existence of similar ha- 
los around external galaxies has been inferred from obser- 
vations of absorption line systems around bright galaxies 
(e.g. Lanzetta et al. 1995). These halos may have suffi- 
cient density to strip any low column density gas mov- 
ing through them. Several investigators have suggested 
that such stripping is responsible for removing gas from 
the Magellanic Clouds as they orbit through the Galaxy's 
halo, producing the purely gaseous Magellanic Stream 
(e.g. Meurer, Bicknell & Gingold 1985; Sofue 1994; Moore 
& Davis 1994). Sofue & Wakamatsu (1993) and Sofue 
(1994) specifically stress that stripping by galaxy halos 
should also play an important role in the evolution of H I 
tidal tails. 

The tidal features in each of our systems have H I col- 
umn densities and velocities similar to those assumed in 
the numerical models of Sofue (1994) and Moore & Davis 
(1994), which resulted in rather extreme stripping of the 
Hi clouds. Although this may seem to provide an im- 
mediate explanation for our observations of gas/star dis- 



placements, we point out that these column densities and 
velocities are typical of all of the tails thus far imaged in 
H I, the great majority of which do not show the extreme 
displacements we describe here. There is no reason to be- 
lieve that the halo properties of NGC 520, Arp 299 and 
Arp 220 are any different from, or that the encounters 
were any more violent than similar mergers which do not 
exhibit such dramatic displacements (e.g., NGC 4038/9, 
NGC 7252, NGC 4676 HvG96; NGC 3628 Dahlem et al. 
1996; NGC 2623, NGC 1614, Mrk 273 Hibbard & Yun 
in preparation; NGC 3256 English et al. 1999). In fact, 
in light of the stripping simulations mentioned above, one 
wonders why such displacements are not more common. 

A possible solution to this puzzle is suggested by the re- 
sults of numerical simulations of major mergers. In these 
simulations, the material distributed throughout the ha- 
los of the progenitors is tidally distended along with the 
tails, forming a broad sheath around them (see e.g. the 
video accompanying Barnes 1992). This sheath has simi- 
lar kinematics as the colder tail material, resulting in much 
lower relative velocities than if the tail was moving through 
a static halo, thereby greatly reducing any relative ram 
pressure force. 

In summary, while halo stripping might be effective 
for discrete systems moving through a static halo (such 
as the LMC/SMC through the halo of the Galaxy, or 
disks through a hot cluster IGM), the lack of widespread 
Hi/optical decoupling in mergers suggests that it is not 
very effective for removing gas from tidal tails, and it does 
not appear to be a suitable explanation of the present ob- 
servations. 

3.4.2. RPS from Expanding Superwind 

The three systems under discussion host massive nuclear 
starbursts with associated powerful outflows or "super- 
winds". Optical emission lines and/or X-ray emission re- 
veal that the observed outflows extend for tens of kpc from 
the nuclear regions. Theoretical calculations suggest that 
the observed gas plumes represent just the hottest, densest 
regions of a much more extensive, lower density medium 
(Wang 1995). In each of these three systems, the most 
extreme gaps in the H I distribution appear along the in- 
ferred direction of the expanding hot superwind (Figs. [j]d, 
[2|d, ^d). A very similar anti-correlation between an out- 
flowing wind and tidally disrupted H I has been observed 
in the M82 system (Yun et al. 1993; Strickland et al. 1997; 
Shopbell & Bland-Hawthorn 1998), the NGC 4631 sys- 
tem (Weliachew et al. 1978; Donahue et al. 1995; Wang 
et al. 1995; Vogler & Pietsch 1996), and possibly the NGC 
3073/9 system (Irwin et al. 1987; Filippcnko & Sargent 
1992). It has been suggested that this anti-correlation is 
due to an interaction between the blown-out gas of the 
superwind and the cold gaseous tidal debris, either as the 
wind expands outward into the debris, or as the tidal de- 
bris passes through the wind (Chevalier & Clegg 1985; 
Heckman, Armus & Miley 1987, and references above). 

Figure || presents the suggested geometry for the case of 
Arp 299. This figure is constructed from our preliminary 
efforts to model the northern tail of Arp 299 using N-body 
simulations similar to those presented in Hibbard & Mihos 
(1995; i.e. no hydrodynamical effects are included). We 
found that we could not match the morphology and kine- 
matics of both filaments simultaneously, but could match 
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Fig. 5. — An illustration of the proposed relative placement of the tidal tail and the expanding wind or ionization cone for the 
case of Arp 299. In this figure, the purely gaseous tidal filament is represented by light grey circles, whereas the optical tidal 
filament is represented by circles drawn with two shades: regions of the optical tail which have accompanying H I are shaded dark 
grey, while t he reg ions which have been cleared of H I are shaded white (see Fig. ^). The suggested geometry of the superwind 
bubble (cf. §|X|) or ionization cone (cf. § p.5.l[ ) is indicated by the cone, which has been drawn to encompass the white circles 
in the optical filament. This figure illustrates how the restricted opening angle of such a cone may only intersect a portion of the 
ribbon-like tail. 



either one separately. Fig. [| presents the results of com- 
bining these two solutions. In this sense this figure is not a 
self- consistent fit to the data, but simply a cartoon which 
illustrates the proposed relative placement of the tidal tail 
and the expanding wind. In this figure, the wind opening 
angle is illustrated by the cone, the gas-rich regions of the 
tails are represented by dark and light grey circles, and 
the gas-poor regions of the tails are represented by white 
circles. The figure illustrates how the restricted o pening 
angle of such a wind (or of an ionization cone, cf. §3.5.1) 
may intersect only a portion of the ribbon-like tail. 

Here we estimate whether ram pressure stripping by the 
nuclear superwind can exert sufficient pressure on gas at 
the large distances typical of tidal tails. We use equa- 
tion (5) from Heckman, Lehnert & Armus (1993; see also 
Chevalier & Clegg 1985) to calculate the expected ram 
pressure (Prps) of the superwind far from the starburst 
as a function of its bolometric luminosity (if, c ;): 



Pnps(r) = 4 x 10 10 dynecm" 



Lboi 



1O"L 



1 kpc 



This equation has been shown to fit the pressure pro- 
file derived from X-ray and optical emission line data of 
Arp 299 (Heckman et al. 1999). It should provide a lower 



limit to the ram pressure, since it assumes that the wind 
expands spherically, while observations suggest that the 
winds are limited in solid angle. 

This pressure can be compared with the pressure of the 
ambient medium in the tidal tail (Pudai), given by the 
energy per unit volume: 



Ptidal — C X PgasO~^ as 



(2) 



where C is a constant, p gas is the mass density of the gas, 
and o-g as is the velocity dispersion of the gas. For an equa- 
tion of state of the form P ~ p 1 , the constant C is equal 
to 7 _1 = |, and a gas corresponds to gas sound speed, 
while for a self-gravitating cloud, C = | and a gas is the 
one-dimensional velocity dispersion of the cloud. In both 
cases, we assume that the observed line-of-sight velocity 
dispersion of the H I is a suitable measure of o~ gas . 

The mass density of the tidal gas is given by p gas = 
1.36 x tuh x tihi, where the numerical constant accounts 
for the presence of He, mn is the mass of a hydro- 
gen atom, and nui is the number density of atomic hy- 
drogen. If we assume the gas is uniformly distributed 6 
with a column density Njji along a length dL, we have 

r )(!2£s), where the 



tihi 



6.5 x 10" 3 — " 3 ' 



Nhi 



1 ^2xl0 20 cm- 2 A dL 

fiducial values are typical of tidal features (HvG96, Hib- 



°Clearly the results will be very different if the tidal gas is mainly in dense clouds, a point that can be tested with higher-resolution VLA 
observations. For now we calculate the ram pressure effect on the diffuse gas. 
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bard & Yun in preparation). We rewrite eqn. 
of the observables: 



(|2|) in terms 



°udai =2x10 14 dyne cm 



C_ 

L5 



N 



HI 



2 x 10 20 cm- 



lOkpc 
dL 



o~hi 



lOkms 



(3) 



The maximum radius out to which we expect material 
to be stripped {Rrps) is then given by the requirement 
that P RP s(r) = Pudali*) at r= Rrps- We replace L bot in 
eqn. (pi) by the IR luminosity (Lip), under the assumption 
that the IR luminosity arises from reprocessed UV photons 
from the starburst 7 (Lonsdale, Persson & Matthews 1984, 
Joseph & Wright 1985). The very high IR luminosities of 
these systems (Ljr/Lb > 10) make it likely that this is 
indeed the case, and we are probably making an error of 
<; 10% (e.g. Heckman, Lehnert & Armus 1993). Equating 
eqns. (P and (||) we find: 



Rrps = 140 kpc 



1.5 
~C 



^IR 



10 n L 



2 x 10 



20, 



N, 



HI 



dL 
10 kpc 



^10kms~ 



<^h i 



(4) 



In Table [l] we provide estimates of Rrps for the sys- 
tems considered here. In calculating Rrps we made the 
very conservative assumption that the values of Nhi and 
(Thi for the stripped gas are equal to the maximum val- 
ues found within the tidal tails (see Table [l]) . The results 
of these calculations indicate that, in all cases, Rrps is 
larger than the radii of the observed gaps in the tidal H I 
distributions. Therefore, in principle, the wind should be 
able to strip the gas from any tidal material in its path. 

The above derivation assumes that the tidal gas is at rest 
with respect to the wind. It can be easily generalized to 
the case of a wind impacting an expanding tidal feature by 
reducing the wind ram pressure by a factor ( ^"^ ) 2 . 
For NGC 520 and Arp 220, the tidal gas is primarily in 
rotation (i.e., moves perpendicular to superwind), so we 
expect the gas to feel the full ram pressure given above. 
For Arp 299, Heckman et al. (1999) find V wtnd = 800 km 
s , and we estimate a maximum Vudal = 240 km s" 1 
(HY99). Therefore, the ram pressure could be reduced b 
50%, reducing Rrps by 70% from that listed in Table 
i.e. Rrps ~ 100 kpc for Arp 299. This is still large enough 
to reach to the region of the N clump in Fig. ||. The regions 
of high H I velocity dispersion indicated in Fig. |^c may be 
due to the influence of such a wind. We note that these 
regions occur on the side of the H I features that face the 
starburst region. However, no such kinematic signatures 
are visible in the gas near the wind axis in NGC 520 or 
Arp 220. 

The lack of gaseous/stellar displacements in the tidal 
tails of many superwind systems might seem to provide 
a strong argument against the scenario outlined above. 
However, there are two conditions needed to produce wind- 
displaced tidal features: the starburst must be of sufficient 



energy and duration to achieve "blowout" , and the tidal 
H I must intersect the path of the expanding wind mate- 
rial. The second condition is not met for the blowout sys- 
tems NGC 4676, NGC 3628, NG C 26 23, NGC 1614, and 
NGC 3256 (references given in §3.4.1). In these systems 
the tidal tails appear to lie at large angles with respect 
to the blowout axis, and their tidal tails should not inter- 
sect the wind. Both conditions are met for M82 and NGC 
4631, which both show extreme H i/optical displacements. 
For these two systems, the high-latitude H I appears to 
be accreted from nearby disturbed companions (M81 and 
NGC 4656, respectively; Yun et al. 1993, Weliachew et al. 
1978), while for the three major mergers under study here 
the H I appears to intercept the path of the wind as a 
result of a highly inclined encounter geometry. A high 
inclination encounter geometry is therefore a prerequisite 
for such displaced morphologies, and the host of the su- 
perwind should be the disk with a retrograde or polar spin 
geometry. 

Nevertheless, ram pressure stripping cannot provide a 
complete explanation of the observations. Since the stars 
are unaffected by the wind, it would be an unusual coin- 
cidence for the edges of the optical plumes to correspond 
with the edge of the cold gas which is presently being ab- 
lated. Nor does it seem likely that the wind could be suf- 
ficiently collimated to "bore" into the northern H I clump 
in Arp 299 just where the optical tail appears projected 
upon it. Therefore a second process is still needed to ex- 
plain the small-scale anti-correlations. In conclusion, the 
expanding wind should affect any tidal H I in its path; how- 
ever this effect alone cannot explain all the details of the 
observations. 



3.5. Photoionization 

Disk galaxies are known to exhibit a precipitous drop in 
neutral hydrogen column density beyond column densities 
of a few times 10 19 cm -2 (Corbelli, Schneider & Salpeter 
1989; van Gorkom 1993; Hoffman et al. 1993). This drop 
has been attributed to a rapid change in the ionization 
fraction of the gas due to influence of the intergalactic UV 
field (Maloney 1993, Corbelli & Salpeter 1993, Dove & 
Shull 1994; see also Felten & Bergeron 1969, Hill 1974), 
rather than a change in the total column density of H. 

Tidal tails are assembled from the outermost regions 
of disk galaxies. Since this material is redistributed over 
a much larger area than it formerly occupied, its surface 
brightness must decrease accordingly. Therefore, if the 
progenitors were typical spirals, with H I disks extending 
to column densities of a few times 10 19 cm" 2 , then the 
resulting tidal tail must have gas at much lower column 
densities. However, tidal tails exhibit a similar edge in 
their column density distribution at a similar column den- 
sity (HvG96; Hibbard & Yun in preparation). This is one 
of the most compelling pieces of evidence for an abrupt 
change in the phase of the gas at low Hi column den- 
sity. The outer tails mapped in H I should therefore be 
the proverbial "tip of the iceberg" of a lower column den- 
sity, mostly ionized medium. With the tidal gas in this 
very diffuse state, fluctuations in the incident ionizing flux 

7 This assumes that the IR, luminosity is not enhanced due to the presence of AGN. There is no evidence for an energetically important AGN 
in any of these three systems. 
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might be expected to produce accompanying fluctuations 
in the neutral gas fraction. 

Given these considerations, we examine the possibility 
that the total hydrogen column density does not change 
at the regions illustrated in Figs. |]-^|, but that the neutral 
fraction does, i.e. that the gas in the regions under study 
has a higher ionization fraction than adjacent regions. The 
intergalactic UV field should be isotropic, and would not 
selectively ionize certain regions of the tails. Here we ex- 
amine the possibility of two non-isotropic sources of ion- 
izing flux: (1) leakage of UV flux from the circumnuclear 
starburst; (2) ionization by late B stars and white dwarfs 
associated with the evolved stellar tidal population. 

Our procedure is to compare the expected ionizing flux 
density shortward of 912 A to the expected surface re- 
combination rate of the gas. We assume that in the 
area of interest the gas is at a temperature of ~ 10 4 
K, for which a case-B recombination coefficient of = 
2.6 x 10~ 13 cm 3 s 1 is appropriate (Spitzer 1956). We as- 
sume that the hydrogen is almost completely ionized, so 
n e ~ hh ■ We further assume that the density of ionized 
gas is the same as the density of the neutral gas in the 
adjacent regions, uh ~ "m, w here uhi is calculated as 
above (n# ~ Njji/dL, § |3.4.2 ). The detailed ionization 
state will depend sensitively on the dumpiness of the gas, 
but a full treatment of this problem is beyond the scope of 
this paper. Here we wish to investigate if these processes 
are in principle able to create effects similar to those ob- 
served. 

3.5.1. Photoionization by the Starburst 

Here we consider the case that the superwind does not 
affect the tidal gas by a direct interaction, but influences 
it by providing a direct path from the tidal regions to 
the starburst, free of dust and dense gas (see also Fig. ||). 
Through these holes, ionizing photons from the young hot 
stars stream out of the nuclear regions and are quickly 
absorbed by the first neutral atoms they encounter. Fol- 
lowing Felton & Bergeron (1969; see also Mahoney 1993), 
we solve the equation: 



O.1M 
yields 9 



and M; 
I 
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fesc (S^"ti 



125Af , and solar metallicity. This 

1.83x10 s4 photonss -1 „ L 
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Again making the standard assumption that most of the 
starburst luminosi ty is emitted in the far infrared (i.e. 
Lboi = L IR , cf. §3.4.2), we rearrange eqn. (|) to solve 
for the radius, Rionized, out to which the starburst is ex- 
pected to ionize a given column density of H I of thickness 
dL: 



R 



ionized 



— 66 kpc 
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dL 
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(0) 



Resulting values for Ri on i ze d are listed in Table [D. For 
this computation, we have adopted a value of 10% for 
fesc(^wind) ■ This is equal to the total fraction of ioniz- 
ing photons, / esc , escaping from a normal disk galaxy as 
calculated by Dove, Shull & Ferrara (1999). Even higher 
values of f esc are expected in starburst systems (Dove et al. 
1999). Since we stipulate that a higher fraction of ionizing 
photons escape along sightlines above the blowout regions 
than are emitted along other directions, it follows that 
fesd^wind) > fesc, and as a result the values of R lon ized 
calculated in Table [l] should be conservative estimates. 

Table |l| shows that under these simplified conditions, 
Rionized is of the order of, or larger than, the tidal radii of 
interest. We therefore conclude that the starburst seems 
quite capable of ionizing tidal H I, if indeed there is an 
unobstructed path from the starburst to the tidal regions. 
This might explain the lack of H I along the wind axis in 
NGC 520 and Arp 220, and the absence of H I along the 
optical tidal tail in Arp 299. 

This process is especially attractive since it can poten- 
tially explain the lack of H I at the bases of otherwise gas- 
rich tidal tails in NGC 7252 (Hibbard et al. 1994), Arp 105 
(Due et al. 1997), and NGC 4039 (Hibbard, van der Hulst 
& Barnes in preparation). These systems do not show 
evidence for expanding superwinds, which rules out the 
possibility that RPS is playing a role. And each of these 
systems possesses a level of star formation that, according 
to eqn. o, is capable of ionizing gas out to the necessary 
radii. 

However, photoionization by the central starburst does 
not seem capable of explaining all of the observations. As 
with the wind hypothesis above, it would be an unusual 
coincidence for the edges of the optical plumes to corre- 
spond with the edge of ionization cone. Therefore a sec- 
ond process is still needed to explain the small-scale anti- 
correlations. 

3.5.2. Photoionization by the Optical Tails 

The fact that the H I column density falls off just as the 
optical surface brightness increases at the edges of various 
tidal features (Fig. [|) leads us to suspect that there may be 
local sources of ionization within the stellar features them- 
selves. For NGC 520 and Arp 220, we believe the outer H I 

8 It is important to differentiate fesc, the total fraction of ionized photons emerging from a starburst, and fesc(^wind), the fraction emerging 
along a particular sightline. f esc is the total angle averaged fraction, i.e. f esc is the integral of f esc (dQ) over all solid angles, while fesc(Q w ind) 
is the integral over a solid angle cleared by the wind. Most studies in the literature quote values for f eS c- 

9 With similar assumptions, the "Starburst99" models of Leitherer et al. (1999) give approximately the same numerical coefficient. 
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= 13.6eV 



The right hand side of this equation represents the to- 
tal ionizing radiation escaping the starburst region along 
a direction that has been cleared of obscuring material 
by the superwind. We express this in terms of the to- 
tal ionizing flux of a completely unobscured starburst of a 
given bolometric luminosity L^oi by introducing the factor 
fesc(^wind) to account for the fact that only a fraction of 
the photons emitted into a solid angle VL W ind nn d their way 
out of the starburst region 8 . The expected ionizing flux for 
a starburst of a given bolometric luminosity Lboi is calcu- 
lated from the population synthesis models of Bruzual & 
Chariot (1993; 1995), assuming continuous star formation 
with a duration longer than 10 Myr (long enough for the 
burst to achieve blowout), a Salpeter IMF with M\ owcr = 



Hibbard, Vacca, & Yun 



13 



is in a disk structure which is intersected by the tidal stel- 
lar plumes and we wish to investigate whether ionization 
by evolved sources within the stellar plumes, such as late 
B stars and white dwarfs, could be responsible for decreas- 
ing the neutral fraction of the diffuse outer H I. For Arp 
299 the geometry is more complicated, and we refer the 
reader to Fig. |^. Here we suggest that part of the purely 
gaseous tidal filament (the light grey filament in Fig. ||) is 
ionized by evolved sources near the end of the stellar tail 
(the dark grey filament in Fig. |5|, especially those regions 
nearest the gas-rich filament in the right hand panel of this 
figure). 

As in the previous section, we balance the surface re- 
combination rate with the expected ionizing flux density 
(eqn. |5|). In this case, we calculate the ionizing flux den- 
sity for an evolved population of stars of a given i?-band 
luminosity density in L Q pc -2 ). 

In order to approximate the stellar populations in the 
tidal tails, we assume that the tails arise from the outer 
edges of an Sbc progenitor, and that star formation ceased 
shortly after the tails were launched. We again use the 
models of Bruzual & Chariot (1993, 1995) for a Salpeter 
IMF over the mass range 0.1-125 M Q , and adopt an ex- 
ponentially decreasing SFR with a time constant of 4 
Gyr (typical of an Sbc galaxy, Bruzual & Chariot 1993), 
which is truncated after 10 Gyr and allowed to age an- 
other 500 Myr. This simulates the situation in which 
star formation within the disk is extinguished as the tail 
forms, and the ejected stellar population passively fades 
thereafter. While tidal tails frequently exhibit in-situ 
star formation (e.g. Schweizer 1978; Mirabel, Lutz & 
Maza 1991), it is usually not widespread. Under these as- 
sumptions, a population with a projected i?-band surface 
brightness of 1 £_r.© pc -2 should produce an ionizing flux 
of 2.36 x 10 4 ph s -1 cm -2 . Therefore eqn. (||) becomes 
n 2 HI aB<iL < 2.36 x 10 4 phs -1 cm 2 x T,r, which can be 
rewritten as: 



E R > 14.28 L fli0 pc- 



Nhi 



2 x 10 20 cm" 2 



lOkpc 
dL 



(7) 



Noting that 1 L Q pc -2 corresponds to ^^=25.9 mag 
arcsec -2 , we rewrite this as a condition on the surface 
brightness of the tidal features: 



fiR < 23.0 mag arcsec 
2.5 x log 



N, 



HI 



2 x 10 20 < 



/ 10 kpc 
V dL 



(8) 



Referring to Fig. |], we see that only the northern plume 
of Arp 220 is bright enough to ionize nearby tidal H I at the 
appropriate column densities. Neither the optical plume 
in the NGC 520 system nor the northern tail in the Arp 
299 system appears bright enough to ionize the necessary 
columns of hydrogen unless the tidal features are unrea- 
sonably thick (~60 kpc). However, since we have no other 
explanation for the small scale Hi/optical differences il- 
lustrated in Fig. ^, we are hesitant to abandon this expla- 
nation too quickly. 

A possible solution is to invoke continued star forma- 
tion even after the tails are ejected. For instance, if we 



do not truncate the star formation rate after 10 Gyr, in- 
stead allowing the star formation rate to continue its ex- 
ponential decline as the tail expands, then the ionizing 
flux per 1 Lr <& pc -2 is 70 times higher than the value 
of 2.36 x 10 4 ph s -1 cm -2 used above. This would lower 
the fiducial surface brightness in eqn. (||) from 23.0 mag 
arcsec -2 to 27.5 mag arcsec -2 , in which case the faint 
tidal features in NGC 520 (fiR ~ 25 mag arcsec -2 ) and 
Arp 299 (fiR ~ 26.5 mag arcsec -2 ) could indeed ionize the 
necessary column densities of adjoining H I. 

The observed broad-band colors of the tidal tails are not 
of sufficient quality to discriminate between these two star 
formation histories, since the expected color differences are 
only of order B — R=0.l mag. However, whether or not 
the gas is more highly ionized in the regions of interest 
can be addressed observationally. The expected emission 
measure (EM — J nidi) can be parameterized as: 



EM = 0.42 cm -b pc 



N, 



HI 



2 x 10 20 cm" 2 



10 kpc 
dL 



(9) 



Since emission measures of order 0.2 cm -6 pc have been 
detected with modern CCD detectors (e.g. Donahue et al. 
1995; Hoopes, Walterbos & Rand 1999), there is some 
hope of being able to observationally determine if regions 
of the tidal tails are significantly ionized. If the gas has 
a clumpy distribution, then there should be some high 
density peaks which might be sufficiently bright to yield 
reliable emission line ratios. Such ratios would allow one 
to determine the nature of the ionizing source, e.g. pho- 
toionization vs. shocks. Therefore, while we cannot assert 
unequivocally that photoionization plays a role in shap- 
ing the outer tidal morphologies, it is possible to test this 
hypothesis with future observations. 

The hypothesis that ionizing flux from a stellar tidal fea- 
ture may ionize gas in a nearby gaseous tidal feature is not 
necessarily at odds with the observations that many stel- 
lar tails are gas-rich. This is because tails with cospatial 
gas and stars arise from regions originally located within 
the stellar disk of the progenitors, while the optical faint 
gas-rich tidal feat ures likely arise from regions beyond the 
optical disk (§ |3.l| ). In normal disk galaxies, the H I within 
the optical disk is dominated by a cooler component with 
a smaller scale height and velocity dispersion, while the 
H I beyond the optical disk is warmer and more diffuse, 
with a larger scale height and velocity dispersion (Braun 
1995, 1997). As a result, dL should be considerably larger 
for purely gaseous tidal features than for optically bright 
tidal features. 

4. CONCLUSIONS 

In this paper we have described differences between 
gaseous and stellar tidal features. There are large-scale 
differences, such as extensive purely gaseous tidal features 
(the outer disks in NGC 520 and Arp 220 and the outer 
filament in Arp 299) and largely gas-poor optical features 
(tidal plumes in NGC 520 and Arp 220 and the inner fila- 
ment in Arp 299). And there are smaller-scale differences: 
the anti-correlation between the edges of gaseous and op- 
tical features depicted in Fig. ||. A similar anti-correlation 
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is observed between Hi and optical shells in shell galax- 
ies (Schiminovich et al. 1994a, b, 1999), many of which are 
believe to be more evolved merger remnants. 

We have examined a number of possible explanations 
for these observations, including dust obscuration, differ- 
ences in the original distribution of gas and starlight in 
the progenitor disks, gas cloud collisions within the de- 
veloping tails, ram pressure stripping due to an extensive 
hot halo or an expanding superwind, and photoionization 
by either the central starburst or evolved sources in the 
tidal tails themselves. However, no one model easily and 
completely explains the observations, and it is conceivable 
that all explanations are playing a role at some level. 

The most likely explanation for the lack of starlight as- 
sociated with the outer tidal H I is that such features arise 
from the H i-rich but optically faint outer radii of the pro- 
genitor disks. The relatively short lifetimes of luminous 
stars and the large velocity dispersions of less luminous 
stars, especially with respect to the gas, will further dilute 
the luminous content of this material, and the H i-to-light 
ratio of the resulting tidal features will increase with time 
(Hibbard et al. 1994). Gaseous tidal extensions with very 
little detectable starlight would seem to be the natural con- 
sequence. The outer H I rings in NGC 520 and Arp 220 
and the gas-rich outer filament in Arp 299 are all likely to 
arise from these gas- rich regions of their progenitor disks. 

For the gas-poor tidal features we suggest that the star- 
burst has played an important role in shaping the gaseous 
morphology, either by sweeping the features clear of gas via 
a high-pressure expanding superwind, or by excavating a 
clear sightline towards the starburst and allowing ionizing 
photons to penetrate the tidal regions. The primary sup- 
porting evidence for this conclusion is rather circumstan- 
tial: the five galaxies with the most striking Hi/optical 
displacements (the three systems currently under study 
here, and the H I accreting starburst systems M82 and 
NGC 4631) host massive nuclear starbursts with associ- 
ated powerful outflows or superwinds aligned with the di- 
rection of the most extreme Hi/optical displacements. 

NGC 520, Arp 299, and Arp 220 each experienced pro- 
grade/polar or prograde/retrograde encounters. This rela- 
tive geometry may be a pre-requisite for the morphological 
differences reported here. Retrograde and polar encoun- 
ters do not raise extensive tidal tails (e.g. Barnes 1988), 
leaving large gaseous disks in the inner regions. These 
disks should help collimate and "mass-load" the superwind 
(Heckman, Lehnert & Armus 1993; Suchkov et al. 1996), 
which in turn leads to denser and longer-lived winds. Si- 
multaneously, the combination of opposite spin geometries 
provides the opportunity for the tidal tail from the pro- 
grade system to rise above the starburst region in the polar 
or retrograde system, where it may intersect the escaping 
superwind or UV radiation. If this suggestion is correct, 
only systems hosting a galactic superwind and experienc- 
ing a high-inclination encounter geometry should exhibit 
such extreme differences between their H I and optical tidal 
morphologies. 

The observations do not allow us to discriminate be- 
tween either the RPS or the photoionization models: sim- 



ple calculations suggest that either is capable of affecting 
the diffuse outer gas if the geometry is right. There might 
be some evidence for the effects of an impinging wind on 
the outer material in Arp 299 from the increased veloc- 
ity dispersion at several points (HY99); however NGC 520 
and Arp 220 show no such signatures. Photoionization 
is an attractive solution, as it offers a means of explain- 
ing the lack of tidal H I found at the base of otherwise 
gas-rich tidal tails in mergers which show no evidence of 
a supe rwind (e.g. NGC 7252, Arp 105, NGC 4039; see 

Since any ionized hydrogen will emit recombination 
lines, both explanations can be checked observationally. 
The expected emission measure is given by eqn. (^]), which 
predicts detectable features at the column densities of in- 
terest. The morphology of the ionized gas should reveal 
the nature of the ionizing source: photoionized gas should 
be smoothly distributed, while gas excited by RPS should 
be concentrated in dense shocked regions on the edges of 
the H I that are being compressed by the superwind, i. e. 
on the edges nearest the wind axis in Figs. [i]d, ^d & ||d. If 
the gas is clumpy, there may be regions bright enough to 
allow line ratios to be measured, which should further aid 
in discriminating between photoionization or shock exci- 
tation. 

Only two scenarios are offered to explain the small-scale 
anti-correlations: dust obscuration and photoionization 
due to evolved sources in the optical tails. Dust obscu- 
ration likely affects the apparent tidal morphologies at the 
lowest light levels, but we suspect that the dust content is 
too low to significantly obscure the brighter tidal features. 
However, if the tidal tails are highly ionized, with the neu- 
tral gas representing only a small fraction of the total hy- 
drogen column density, it is possible that we are grossly 
underestimating the expected amount of absorption. This 
question can be investigated directly with submm imaging 
of the cold dust in tidal tails. 

The other possibility is that the UV flux from evolved 
sources in the optical tails is responsible for ionizing nearby 
diffuse outer Hi. A simple calculation suggest that the 
tidal tail in Arp 220 is bright enough to ionize nearby 
H I, but the expected ionization flux from the optical tails 
in NGC 520 and Arp 299 is too low to explain the ob- 
served differences, unless significant star formation con- 
tinued within these features after their tidal ejection. If 
this is indeed the case, then the regions where the neutral 
gas column density drops rapidly (see Fig.||) should con- 
tain ionized gas which would emit recombination radia- 
tion. The expected levels of emission should be observable 
with deep imaging techniques (see above). This situation 
requires that the gas and stellar features are physically 
close, and not just close in projection, which can tested 
with detailed numerical simulations. 
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Table 1 

Effects of Ram Pressure and Ionization on Tidal Material 



Parameter Units NGC 520 Arp 299 Arp 220 notes 



V he i (kms- 1 ) 2260 3080 5400 

D (Mpc) 30 48 79 a 

L IR (L ) 7.6 x 10 10 8.1 x 10 11 1.5 x 10 12 b 

N m , m ax (cm- 2 ) 5 x 10 20 4 x 10 20 4 x 10 20 c 

(THi'max (kms^ 1 ) 16 20 40 c 

r (kpc) 25 70 30 d 

R RPS (kpc) 50 140 95 e 

Rionize (kpc) 25 95 130 / 



a Adopting the distances from Sanders, Scoville & Soifer (1991), which assumes H Q =75 km s 1 Mpc 1 and the Virgo- 
centric flow model of Aaronson et al. (1982). 

b 8 — 1000 fim luminosity, from Sanders, Scoville, & Soifer (1991). 

c Maximum values of gas column density and line-of-sight velocity dispersion observed within the tidal features, taken 
from HvG96 (NGC 520), HY99 (Arp 299) and Yun & Hibbard 1999a (Arp 220). 

d The projected radius of "missing" H I. 

°Calculated via eqn. (|]), assuming the values given in this table and <iL=10 kpc. 

f Calculated via eqn. (^j), assuming the values given in this table and f e sc{^wind)—0.1 and c£L=10 kpc. 



